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Thermal  Loads  Are  Rising  Sharply 


•  Thermal  load  in  satellites  is  doubling  in  every  5.5  years 

•  Space  Structures 

-  Heat  flux  of  700  W/cm2  (hot  spots  by  actuators,  etc.) 

•  To  spread  the  hot  spots  (AT~1°C  across  Ax~1cm) 
requires  material  of  thermal  conductivity  of  k~70 
W/m-K 

•  In  comparison  k  for  adhesive  -  0.3  W/m-K 


•  Electronics  Cooling 

-  ~  60-70  W/m-K  is  desirable  for  Electronics  Heat 

Sink  system 


Heat  exchanger 


•  Need  for  high  fidelity  thermal  component  design  for 
tailoring  its  thermal  properties  to  meet  system 
requirements 


Electronic  cooling 


Heat  Sink  Interface 
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Management  of  Thermal  Energy  in 
Materials  &  Systems 


Thermal  Energy  Mgmt 


Thermal 
transport  in 
materials  and 
system 
components 

Thermal  energy 
storage 

Thermal  energy 
conversion 

Etc. 


Technical  Approaches 

•  Passive  system 
(tailoring 
material  thermal 
properties) 

Active  system 
icro  porous 
heat  fluid  flow, 
etc.) 

Etc. 


Thermal  materials  & 
its  interface 
property  tailoring 


Technical  Challenge 


•  Numerous  prior  efforts  have  been  made  by  mixing 
of  CNT  in  polymers  (epoxy)  yields  limited 
improvement  in  thermal  conductivity  (k) 

•  k  (SWNT)  -  2000  -  6000  W/m-K 

•  Improvement  is  limited  to  only  125%  (k  -  0.7  W/m-K) 
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V 


JAccV  Spot  Magn  Det  WD  Exp  I — 

|l00kV  30  5000X  SE  7  5  2  I M7/PR24  PS/862/4%  CROSS  SEC 


"T  5  pm 


M.  J.  Biercuk,  et  al,  Apply.  Phys.  Lett.,  80,  2767  (2002) 


Primary  reason  of  the  limited  improvement 
*  Phonon  scattering  at  the  CNT-polymer  interface 


Technical  Challenge 

*  Thermal  interface  design  for  aerospace  materials 

—  Nano  constituents  interface  in  presence  of  amorphous  materials 
(composites  and  adhesive  joints,  etc.) 
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Two  Examples  of  Thermal  Materials 
Design  for  Aerospace  Systems 


Hybrid  Fiber  Composites 


•  Intermingled  network  of  Nano 
platelets  grown  on  carbon 
fibers  -  embedded  in  polymer 


o  200  4QQ  m  m  1000 

Pitch  k_l_NT  f\\7m~K] 

fiber 

Roy,  et  al.,  Diamond  &  Related  Materials,  2009 


Aligned  CNTs  in  adhesive  joints 

•  Incorporation  of  MWNTs  along 
the  thickness  of  the  adhesive 
joint 


Press  &  Cure 


conductive  transition 
zone  between  MWNT 
and  adherents 


Outline 


•  Molecular  dynamics  simulation  of  thermal 
transport  in  cross-linked  polymers 

•  Comparison  of  various  energy  components  in 
polymer  thermal  transport 

•  Thermal  interface  resistance  of  CNT/polymer 
interface 

•  Thermal  property  measurement 

-  Characterization  tools  under  development 

-  EELS  technique 

-  Micro  heater 


Heat  Transport  Modeling  of  Epoxy 

Networks 


Calculation  of  Thermal  Conductivity 


Green-Kubo  Approach  (Equilibrium  MD):  This  approach  uses  concept  of  fluctuation-dissipation  theorem 
which  relates  equilibrium  fluctuations  to  out  of  equilibrium  properties  via  an  autocorrelation  function 
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Fourier  Approach  (Non-Equilibrium  MD):  This  approach ,  also  known  as  direct  method,  is  analogus  to 
experimental  measurement.  It  is  based  on  the  principle  that  heat  flux  at  certain  cross-section  is  directly 
proportional  to  temperature  gradient  at  that  surface. 


dT/dx  =  Temperature  gradient  ^ 

dQ/Adt  =  Heat  flux  per  unit  area  per  unit  time  A  x 


Cold  Hot  Cold 


U -  Lz  - W 
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Heat  flow  per  unit  area  (Kcal/mol-Angs“) 


NEMD  Simulations:  Thermal  Conductivity 


Cold  Region 


Hot  Region 


Cold  Region 


Z  = 


dQ /  dt 


Thermal  conductivity  of  the  cross  linked  network  was  found  to  be  ~0.3  W/m-K  which  is  in  nice 

agreement  with  experimental  findings. 
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Normalized  Heat  Flux  Autocorrelation  Function 


D  Simulations:  Thermal  Conducti 
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2  ^  Comparison  between  both  approaches 
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Experimental  values  of  thermal  conductivity  of 
epoxy  networks  is  ~0.28  W/m-K. 
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Energetic  Contributions  to  Thermal 

Conductivity 
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1:  Contribution  due  to  kinetic  energy; 

2,  3:  Contribution  due  to  potential  energy  (vdwl  and  electrostatic ),  respectively 
4,  5:  Contribution  due  to  short  range  forces  (vdwl  and  electrostatic),  respectively 
6:  Contribution  due  to  long  range  forces  (electrostatic  interactions:  Ewald 

Sum) 

7,  8:  Contribution  due  to  bonded  interactions  (bond  stretching  and  angle  bending) 


Terms  1,  2  and  3  are  known  as  Convective  contributions. 
Rest  of  them  are  known  as  Virial  contributions. 
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Energy  Contribution  Analysis 


Overall  Contributions  Virial  Contributions 

■  Virial  contributions  M  Convective  contributions  ■  vdwl  forces  W  Electrostatic  forces  w  Bonded  forces 


Convective  Contributions 

■  Kinetic  energy  hi  Long  range  electrostatic  energy 

M  Short  range  potential  energy 


Long  Range  Electrostatic  Force 
Contribution 

■  Rest  M  Long  range  electrostatic  forces 


Virial  ( collision )  contribution  is  significantly  larger  than  convective  terms. 

Van  der  Waals  interaction  and  corresponding  forces  are  the  dominant  contributors  for 
thermal  transport  in  polymers. 

Electrostatic  and  bonded  contributions  are  negligible. 
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Interface  Thermal  Conductance  at  CNT 

Epoxy  Interfaces. 


System  Studied 


Heat  Input  at  the  constant  rate  in  to  CNT  in  the 
center 

Heat  Extraction  from  the  outermost  depicted 
(Blue)  shell 


We  have  used  previously  discussed 
algorithm  (with  minor  alterations)  to  build 
shown  nano-composite  system  with 
functionalized  nanotubes. 


Nanotube  functionalized  with 
DETDA 


^  QueatFlux  ^  ''interface 
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Systems  Studied 


^  C^Heat  Flux  per  unit  length  ^  ^  ^Interface 


The  arrow  depicted  the  direction  of  heat  flow. 


Two  types  of  simulations  were  performed. 

a)  Epoxy-freezed 

b)  Epoxy-moving 


1  Link 


Links 


10  Links 


CH2  Linkages 


2  Links  5 


Other  Studied  Cases 


CJ14  Link  CjH6  Link  C4H8  Link 


Cs^Link 


EPON  Link  Aromatic  Link 
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Results 


Temperature  evolution 


Steady  state  temperatures  vs. 
CH2  linkages 


Red:  Temp,  evolution  of  hot  nanotube 
Blue:  Temp,  evolution  of  cold  nanotube. 

The  orientation  of  nanotubes  with  respect  to 
system  dimensions  (top  view)  is  schematically 
shown  in  the  inset. 


Solid  Lines:  Epoxy  Moving  simulations 
Dashed  Lines:  Epoxy  Freezed  simulations 

Inset  shows  the  temperature  drop  between  the  two 
nanotubes  as  the  function  of  CH2  linking. 
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Interface  Thermal  Resistance  across  CNTs: 

Transverse  Connection 


Table  1:  Campari  son  of  Interface  Thermal  Resistance 

I,  It  era  lure  Studies 

Thermal  Interface  Resistance  (m2-K/W) 

Gao  et  ai  ( C  NT- o i  l/po Ly  me  r J 

1.58  x  10“ 8  (ThconrtLcal)40 

Foyg^L  et  aL 

lO^-  lO*  K/W  (Theoir  tlcal)  38 

Brining  et  ai  (CNT-Epoxy) 

0.24-2,6  x  10  8  (Experimental)41 

Huxtable  €t  ai  (CNT-water  suspensions) 

S  3  x  10“ 8  (Experimental)11 

Cola  et  ai. 

0.2-7  x  10" 8  (Experimental) 15 

Huxtable  et  ai  (CNT-octane) 

4.0  x  10“ 8  (Simulation)11 

Shenogin  etai  (CNT-octane) 

3,3  x  I0~8  (Simulation)21 

Clancy  et  ai  (CNT-polymer) 

0,2-9. 6  x  I0-B  (Simulation)22 

Carlborg  er  ai  (CNT-argon) 

40.0,  62.5  x  10  8  (Simulation)25 

Murayania  et  ai 

6,5  x  I0“B  (Simulation)42 

Zhongef  ai 

3,0-12.0  x  10  s  (Simulation) 26 

Xu  tr/  aL 

0.01-0.5  x  I0-B  (Simulation)27 

Current  Study 

0.7- 2, 9  x  10“a  (Simulation) 

Effect  of  linkage  length  as  well  as  their  no. 
on  overall  interface  conductance 


Interface  conductance  increase  with  number  of  linkages  but  decreases  with  overall  linker- 
length. 

The  deviation  from  linearity  at  higher  degree  of  functionality  is  attributed  to  the  inter¬ 
linkages  interactions. 

Two  predicted  values  in  the  table  are  based  on  use  of  two  different  surface  widths  for  cross- 
sectional  areas  consideration,  a)  3.4  Angs  and  b)  Nantube  diameter. 


Wave  Packets  Analysis 
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How  to  measure  the  energy  transmission 


z  [a] 


Approved  for  Public  Release,  Distribution  A:  88ABW-201 0-3277 


Single  mode  phonon  transmission 
in  functionalized  CNT 


1.  Effect  of  the  functionalization  on 
the  phonon  energy  transmission. 


-CH3 

(16  u) 


■C2H3 

(27  u) 


2.  Effect  of  the  difference  in  the 
functional  group .... 


-COOH 

(45  u) 


Longitudinal  Acoustic  mode 
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Dispersion  relation:  (6,6)  SWCNT 


hv  [meV] 
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[km/s] 


Point  Mass  defects: 
Longitudinal  Acoustic  mode 


hv  [meV] 
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[km/s] 


Group  velocity  effect: 
Longitudinal  Acoustic  mode 
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Coupling  to  the  functional  group 


25 

20 

15  ^ 

£ 

10  ^ 


hv  [meV] 


o 

PL, 

3JL0 


10 


OX) 


[km/s] 


CNT  with  vertical -CH 2  linkers 


Cross-section  of  ( 6, 6)  ncmotuhe 


12  linkages 


6  linkages 


3  linkages 


Phonon  Energy  Transmission 
in  CNT  with  vertical  links 


iJL P 


4  Accoustic  polarizations  and  the  radial 
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Electron  Emission  Loss  Spectroscopy 
(Orbital  Picture  of  Ethylene) 


http:/ /www.  uaf.  edu/chem/green 

Chapterl  McMurry.ppt 


EELS  Spectrum  Analysis 


•  pi  -  pi*  transition 

•  sigma  -  sigma*  transition 

•  pi  -  sigma*  interbrand  transition 

•  sigma  -  pi*  interbrand  transition 

•  In  HOPG 

pi  transition  -  ~  6  eV 
sigma  transition  ~  27  eV 
In  diamond 

sigma  transition  ~  34  eV 
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COOH  Nanotubes  in  Epoxy  Matrix 


Nanotu 


Epoxy 


Interface 
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2  nm 


STEM 

Image 


Post  EELS  Analysis 


Counts 


Nanotube  EELS  Spectra 
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Counts 


Epoxy  EELS  Spectra 


sigma 
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Counts 


Pure  Nanotube/Epoxy  Interface  EELS 


Energy  Loss  (eV) 
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Counts 


COOH  Nanotube/Epoxy  Interface  EELS 


sigma 


26.8  eV 


Energy  Loss  (eV) 
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Sub-micron  Scale  Thermal  Conductivity 

Measurement 


Resistance  Temp  Detector  -  RTD 


Another  test 
configuration 
using  the  RTD 


Versatile  RTD  design  for  nano-  to  sub-micron  scale 
direct  thermal  conductivity  measurement 
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Summary 


Thermal  transport  mechanism 
in  amorphous  materials 

£ 

systems  I 

-  Non  bonded  interaction  provides 
the  most  energy  to  thermal 
transport  in  amorphous  materials 
systems 
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Bonds  between  carbon  atoms 
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-  Interface  covalent  bonding 
between  polymers  and  nano 
constituents  surfaces  is  a 
necessity  for  improving  interface  k 

•  Phonon  wave  packet  dynamics  to 
visualize  the  phonon  scattering  &  to 
calculate  the  transmission  function  in 
meso  scale  heat  transfer 
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